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ABSTRACT 

We observed the coronally active eclipsing binary AR Lac with the High Energy Transmission 
Grating on Chandra for a total of 97 ks, spaced over five orbits, at quadratures and conjunc- 
tions. Contemporaneous and simultaneous EUV spectra and photometry were also obtained 
with the Extreme Ultraviolet Explorer. Significant variability in both X-ray and EUV fiuxes 
were observed, dominated by at least one X-ray flare and one EUV fiare. We saw no evidence 
of primary or secondary eclipses, but exposures at these phases were short and intrinsic vari- 
ability compromised detection of any geometric modulation. X-ray fiux modulation was largest 
at high temperature, indicative of fiare heating of coronal plasma rather than changes in emit- 
ting volume or global emission measure. Analysis of spectral line widths interpreted in terms 
of Doppler broadening suggests that both binary stellar components are active. Based on line 
fiuxes obtained from total integrated spectra, we have modeled the emission measure and abun- 
dance distributions. The EUV spectral line fluxes were particularly useful for constraining the 
parameters of the "cool" (< 2 x 10^ K) plasma. A strong maximum was found in the differential 
emission measure, characterized by two apparent peaks at logT = 6.9 and 7.4, together with a 
weak but signiflcant cooler maximum near logT = 6.2, and a moderately strong hot tail from 
logT = 7.6 — 8.2. Coronal abundances have a broad distribution and show no simple correlation 
with first ionization potential. While the resulting model spectrum generally agrees very well 
with the observed spectrum, there are some significant discrepancies, especially among the many 
Fe L-lines. Both the emission measure and abundance distributions are qualitatively similar to 
prior determinations from other X-ray and ultraviolet spectra, indicating some long-term stability 
in the overall coronal structure. 

Subject headings: stars: coronae — stars: abundances — stars: individual (AR Lac) — X-rays: stars 

1. Introduction eclipses can help constrain active region geometry, 

it has been a key system for studying the structure 
AR Lac (HD 210334, HR 8448) is one of the photospheric spots from visible light modula- 

brightest totally eclipsing RS CVn binaries. Since 



tion, the chromosphere from emission of magne- 
sium, calcium and hydrogen, the transition region 
through ultraviolet emission lines, and the coronae 
via emission at extreme ultraviolet (EUV), X-ray, 
and radio wavelengths. There is as yet no compre- 
hensive predictive theory that explains in detail 
all the aspects of coronal emission based only on 
fundamental stellar parameters. Observational at- 
tack is then aimed at increasing the quantity and 
quality of spectral and photometric data to pro- 
vide insights and to help constrain the dependence 
of coronal activity on stellar evolutionary param- 
eters, other magnetic activity indicators or as yet 
unidentified parameters. With the High Energy 
Transmission Grating Spectrometer (HETGS) on 
the Chandra X-Ray Observatory, we are able to 
greatly improve the quality of X-ray spectra by re- 
solving a multitude of coronal emission lines due 
to iron and the hydrogen-like and helium-like lines 
of a number of abundant elements. These spec- 
tra provide both line and continuum fluxes and 
their time variability. The aim of this study is 
to model the lines and continuum with the latest 
atomic data in order to determine the coronal tem- 
perature structure, density and absolute elemental 
abundances. 

AR Lac is comprised of G- and K-type sub- 
giants in a 1.98 day orbit. The components 
are each of slightly greater than one Solar mass 
(1.23M0 and 1.27Mq, respectively) and have radii 
of 1.52i?0 and 2.72Rq. They reach a maximum 
radial velocity separation of 230kms~^ and have 
rotational velocities of 39 and 70 kms~^. At a 
distance of 42 pc, AR Lac is relatively bright. 
Gehren, Ottmann & Reetz (1999) summarized 
these and other fundamental properties of AR Lac. 

AR Lac was detected in X-rays by HEAO 1 
at a luminosity of ^ 10"^"'^ ergs s^-'^ (Walter et al. 
1980); Walter, Gibson & Basri (1983) present an 
early analysis of the X-ray coronae in which it 
was inferred that emission arose from both stellar 
components in compact and extended structures. 
Subsequent X-ray studies were undertaken by 
Ottmann, Schmitt & Kuerster (1993) {ROSAT), 
White et al. (1990) (EXOSAT), Rodono et al. 
(1999) (Beppo-SAX), and White et al. (1994) 
{ASCA). 

Observations with RO SAT and ASCA detected 
a deep primary eclipse and a smaller secondary 
eclipse (Ottmann & Schmitt 1994; White et al. 



1994), while EXOSAT and EUVE observations 
could only confirm the primary eclipse due to 
flares or instrumental limitations (White et al. 
1990; Christian et al. 1996; Brickhouse et al. 1999; 
Pease et al. 2002). Some extended chromospheric 
material was also detected by Montes et al. (1997) 
and Frasca et al. (2000). Analysis of high resolu- 
tion EUVE data (Griffiths & Jordan 1998; Brick- 
house et al. 1999; Sanz-Forcada, Brickhouse & 
Duprce 2003) showed a corona dominated by ma- 
terial at T ~ 10^-^ K, and a substantial amount of 
material at hotter temperatures of T ~ 10^'^ K. 

2. Observations and Data Processing 
2.1. Chandra/H.I]TGS 

We observed AR Lac six times (observation 
identifiers (OID), 6-11) with the Chandra X- 
Ray Observatory's HETG/ACIS-S instrument 
(HETGS) (Weisskopf et al. 2002) in standard 
timed-event mode. Two longer exposures (35 
ks) were at the same orbital quadrature. Four 
shorter observations (~ 8 ks) were paired at the 
two eclipses.^ In an attempt to minimize uncer- 
tainties from any long term trends in activity, 
all the data were taken within five orbits of this 
two day period binary. We applied the ephemeris 
of Ferryman et al. (1997), which differs from the 
more recent determination by Marino et al. (1998) 
by less than 0.01 in phase at the epoch of obser- 
vation. Observational details are given in Table 1. 

The event files were re-processed to apply up- 
dated calibration files (namely, CCD gain and 
bad-pixel filters; we used ASCDSVER CIAO 2.1 
Wednesday, February 28, 2001 and CALDBVER 
2.3). Data were also "de-streaked" to filter out 
the instrumental artifact on CCD-8 and then 
processed to grating coordinates. These event 
lists were then binned onto the standard CIAO^ 
HETGS spectral grids. Effective area tables (aux- 
iliary response files, or ARFs) were made with 
CIAO software (mkgarf ) for each of the HEG and 
MEG gratings, and for -|-1 and —1 orders. We 
show a summary flux spectrum in Figure 1. 

-'Caveat: due to an 0.5 day ephemeris error, comments on 
the phase in the standard data product headers are wrong, 
^http : cxc . harvard . edu 
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2.2. EUVE 



3.2. Variability 



EUVE spectrographs cover the spectral range 
70-180 A, 170-370 A and 300-750 A for the short- 
wavelength (SW), medium- wavelength (MW) and 
long-wavelength (LW) spectrometers respectively, 
with corresponding spectral dispersion of AA ~ 
0.067, 0.135, and 0.270 A/pixel, and an effec- 
tive spectral resolution of 200 400. The 
Deep (DS) Survey Imager has a band pass of 80- 
180 A (Haisch, Bowyer & Malina 1993). Stan- 
dard data products from the EUVE observations 
of AR Lac were obtained through the Multimission 
Archive at Space Telescope (MAST), correspond- 
ing to three observational campaigns starting in 
1993 October 12 (96 ks), 1997 July 3 (74 ks) and 
2000 September 14 (63 ks). 

3. Photometric Analysis 
3.1. Lightcurves 

We made light curves of the Chandra data using 
the CI AO program light curve, and filtered the 
input events so that —3 to +3 orders (excluding 
zero, which is piled'"*), and MEG and HEG pho- 
tons within 1 25 A were all binned into one curve 
for each OID. We also made light curves in some 
stronger line and continuum bands by filtering on 
narrow wavelength regions. We show some of the 
light curves in Figure 2. 

EUVE light curves (small open squares in Fig- 
ure 2) were built from the DS image by taking a 
circle centered on the source, and subtracting the 
sky background within an annulus around the cen- 
ter. Standard procedures were used in the IRAF 
package EUV v. 1.9, with a time binning of 600 s. 
In the following discussion and analysis, the obser- 
vations in 1993 and 1997 are referred to as "quies- 
cent" states; despite of the presence of some minor 
flares the apparently quiescent component dom- 
inates the integrated flux (Sanz-Forcada, Brick- 
house & Dupree 2003). The observations in 2000 
are instead referred to as a "flare" state. Only the 
observations of the 2000 campaign, contemporane- 
ous with the Chandra observations, are employed 
here. 



' "Pilcup" refers to coincidence of photons in the same spa- 
tial and temporal bins. It makes the response non-linear, 
and also, if severe, censors counts via on-board filtering. 



From the X-ray light curve in Figure 2 we see 
that AR Lac was highly variable during the obser- 
vations. There is one obvious flare near phase 0.5 
(orbit 2.5; OID 7) when the count rate increased 
by nearly a factor of four, and then rapidly de- 
cayed. A little more than half an orbit later (orbit 
3, phase 0.13-0.33, OID 9), the count rate was 
decreasing. It is tempting to associate this with 
the tail of the larger flare, but it could just as 
well be an independent event. For the two con- 
secutive segments to connect smoothly, the decay 
rate would have to be variable in a complicated 
way; the estimated decay rate of OID 9 is too 
short to connect with the previous flare. Given 
that OID 8 is somewhat elevated in count rate 
and fairly steady, this would be consistent with 
the high flare frequency and broad range of decay 
rates that AR Lac typically exhibits, as revealed 
by the extensive EUVE and Chandra lightcurves 
analysed by Pease et al. (2002). 

We can use the EUVE light curve (Figure 2) 
to provide a context for the disjoint X-ray obser- 
vations. The Chandra OID 7 flare was simultane- 
ously observed with the EUVE, which showed a 
small step up in count rate. A large EUV flare 
occurred between OID 6 and 8. Assuming that 
the high X-ray to EUV enhancement ratio seen in 
OID 7 always holds, which would be more indica- 
tive of further heating of an existing volume of hot 
plasma rather than evaporation of cold chromo- 
spheric material, the EUV flare must have started 
after OID 6. The basehne trend from OID 8-7-9 is 
probably indicative of the decay of the large flare. 
There is no obvious quiescent state; the nearest 
hint of one is from phase 0.25-0.30. One observa- 
tion (OID 9) is decreasing to this level, but the 
other (OID 6) is actually increasing from a lower 
flux state. 

There are no obvious X-ray eclipses, but ex- 
posures were short at these phases. The dotted 
curve in Figure 2 is a simple occultation model 
with uniform disks of equal surface brightness and 
with relative radii in proportion to the AR Lac 
components, scaled arbitrarily. The lower X-ray 
flux state during primary eclipse egress (OID 11; 
phase 0.0-0.1) is constant. This implies that the 
emerging, smaller G-star is X-ray dark, is dark on 
the portion being exposed, or that emission struc- 
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tures are large and polar so they are not occulted. 
We will rule out the former case later, based on 
line-profile information (see 4.4). 

To investigate qualitative temperature changes, 
we binned light curves in narrow bands, for both 
continuum and strong line regions. Continuua 
near lines were used to derive net line rates, and 
some features with similar temperatures of max- 
imum emissivity were summed to improve statis- 
tics. An example is shown in Figure 2 for Si XIV 
(6.18 A, logTmax = 7.4), which follows the over- 
all integrated light curve, and Si xiii (6.6-6.8 A, 
logTmax = 7.0), whose not rate shows little varia- 
tion. As a qualitative temperature diagnostic, we 
computed the modulation in the count rate, r, for 
each feature, defined as (rmax-^-min)/ (^nax+T-min), 
which is 1.0 for maximum modulation {r^i^i — 0), 
to 0.0 for constant count rate (r^ax = ''min)- In 
Figure 3 we show the modulation as a function of 
temperature of formation. For lines, the tempera- 
ture of formation is defined as the temperature of 
maximum emissivity. 

The trend in modulation is clear from these 
light curves: higher temperature emission is more 
modulated. The cutoff is quite sharp for the 
lines at about logT = 7.0. The ion temperature, 
though, is not a unique diagnostic. For exam- 
ple, the difference in modulation between Mg xii 
and Si xiil is not contradictory, since Mg xii is 
hydrogen-like and has a long emissivity tail ex- 
tending to higher temperatures. The continuum 
modulation showed a similar trend, in that the 
continuum flux at shorter wavelengths, which is 
formed by higher temperature plasma, is also more 
strongly modulated. 

As another qualitative diagnostic of tempera- 
ture changes, we looked at band light curve ra- 
tios to sec if the decrease in OID 9 is a flare de- 
cay (after the EUVE observation, on the extreme 
tail of the EUV flare), and whether the slight rise 
in OID 6 in the same phase interval (but before 
the EUV flare) could be due to rotational modula- 
tion of asymmetrically distributed X-ray emitting 
structures. We used the strongly modulated 1.9- 
2.9 A band light curve for the high temperature 
diagnostic, the relatively unmodulated Fe xvii re- 
gions' continuum bands (14.7-14.9, 16.4-16.6, and 
17.2-17.5 A) as well as Ne x to sample the cooler 
plasma. The ratios clearly showed the flares near 
phases and 0.5, but no signiflcant differences at 



the quadrature phases. We cannot say whether 
either the slow rise or fall are due to flares or rota- 
tional modulation. But given the large EUV flare, 
the OID 9 decrease is probably due to flare decay, 
and temperature changes are below the sensitivity 
of the observations. 

4. Spectroscopic Analysis 
4.1. Spectral Line Fluxes 

Chandra Spectra 

Obtaining line fluxes from coronal X-ray spec- 
tra such as that of AR Lac is complicated by the 
presence of continuum flux. This continuum has a 
shape which is dependent primarily on the plasma 
temperature. In order to extract line fluxes, we 
therefore adopted an iterative approach, whereby 
temperature information from spectral lines was 
used to calculate a model continuum, which was 
then used to refine measurements of the line fluxes. 

We summed ±P* orders and flt line fluxes with 
the CXC software suite, ISIS (Houck & Deni- 
cola 2000).^ Emission lines were fit by convolv- 
ing intrinsic source line profiles (Gaussians) plus 
a model plasma continuum by the instrumental 
response or line spread function (LSF). The free 
parameters were the Gaussian centroids and areas 
of each line, and, if necessary, the Gaussian dis- 
persion. Typically, the lines are unresolved, so the 
Gaussian dispersion was frozen at a value well be- 
low the instrumental resolution. For some cases, 
either due to blends or to orbital velocity sepa- 
ration, we flt the line width (line widths are dis- 
cussed in 4.4). The redistribution component of 
the response for grating spectra is the line spread 
function (LSF) , which is stored in grating redistri- 
bution matrix files (RMF) in the CIAO calibration 
database. HEG and MEG spectra were kept sepa- 
rate, but fit simultaneously. The continuum model 
was obtained iteratively, first by fitting relatively 
line-free regions with a single temperature plasma 
model, and subsequently by using the result of dif- 
ferential emission measure models to improve the 
predicted continuum. The continuum normaliza- 
tion was not allowed to vary in the fitting since the 
apparent continuum is often significantly above 
the true continuum due to line blending. For the 

"'http : //space .mit . edu/CXC/ISIS 



4 



continuum model, wc use the summed true con- 
tinuum and pseudo-continuum components in the 
Astrophysical Plasma Emission Database (APED; 
Smith et al. 2001). 

To be explicit, if the model line contribution, 
5, to a region is expressed as a sum of normalized 
Gaussians, g, enumerated by component i, as 



(1) 



then the n*^ iteration predicted counts spectrum 
in the region of interest can be defined as 

C{h)= [ dXR{X,h)A{X)[S{X) + S^{X)]. (2) 

Here, h is the detected channel, A the wavelength, 
Sc the continuum source model, R the redistri- 
bution function (LSF; or response matrix, RMF), 
and A the effective area. The Ui values are the line 
fluxes, and Aj and trj are the line centroid and dis- 
persion, respectively. The line fluxes determined 
by minimizing this against the counts are listed in 
Table 2. 

4.1.2. EUVE Spectra 

Spectra of AR Lac, binned over the total ob- 
servation, were used to provide fluxes of lines in 
the range ^90-400 A as listed in Table 2. To cor- 
rect the observed fluxes for interstellar hydrogen 
and helium continuum absorption, wc used a ratio 
He i/H 1=0.09 (Kimble et al. 1993), and a value 
of Nh — 1.8 X lO^'^, calculated from the Fe XVI 
A335/361 hue ratio. Lines of Fe ix-Fe xxiv (ex- 
cept Fe xvii) are formed in this spectral range, 
providing good coverage in temperature from only 
one element over the range T ~ 10^'^ — 10^ ''' K, 
avoiding the introduction of further uncertainties 
from the calculation of abundances. 

4.2. Temperature Distribution 

A more quantitative description of plasma tem- 
peratures is given by the emission measure distri- 
bution. This is a one-dimensional characteriza- 
tion of an emitting plasma describing the emit- 
ting power as a function of temperature. It does 
not tell us how the material is arranged geomet- 
rically. This must be derived from other informa- 
tion, such as eclipse or rotational modulation, or 
inferred through (hopefully realistic) assumptions 



of parametric or semi-empirical models such as hy- 
drostatic, magnetically confined loops. Nonethe- 
less, the emission measure distribution remains a 
useful quantity for visualization and comparative 
study of coronal temperature structure. Bowyer, 
Drake & Vennes (2000) present a good review of 
emission measure modeling in the context of cool 
stars and extreme ultraviolet spectroscopy. 

There are many pitfalls in inverting the emis- 
sion integral (Craig & Brown 1976; Hubeny & 
Judge 1995; Mcintosh, Brown & Judge 1998; 
Judge, Hubeny & Brown 1997). Inversion is lim- 
ited by the emissivity functions of real atoms, 
which do not form an orthogonal set of basis func- 
tions; some kind of regularization is necessary. 
Kashyap & Drake (1998) also discuss problems 
of spurious structure caused by errors in atomic 
data, and describe an approach to the problem 
based on a monte carlo technique converged by a 
Markov chain which includes atomic data as well 
as line measurement uncertainties. 

In this analysis, we fit the differential emission 
measure (DEM) and abundances simultaneously 
by minimizing the integrated line flux residuals us- 
ing the emissivities from the APED and the ioniza- 
tion balance of Mazzotta et al. (1998). Our basic 
method is described by Huenemoerder, Canizares 
& Schulz (2001). We minimize a statistic. 



E 

1=1 



1 



AlogT 



N 

t=i 

(3) 

Here, I is a spectral feature index, and t is the tem- 
perature index. The measured quantities are the 
line fluxes, fi, with uncertainties ai. The a priori 
given information are the emissivities, e;t, and the 
source distance, d. The minimization provides a 
solution for the differential emission measure, Dt 
and abundances of elements Z, Az- The exponen- 
tiation of \nD is simply a trick which forces Dt 
to be non-negative, and ^' is a smoothing opera- 
tor which imposes some implicit regularization on 
the solution; we use a Gaussian convolution with 
a dispersion of 0.15 dex. We omit spectral fea- 
tures which are line blends of different elements 
and of comparable strengths. We use a temper- 
ature grid of 60 points spaced by 0.05 in logT, 
from logT = 5.5 to 8.5. The emissivities, e, are 
as defined by Raymond & Brickhouse (1996). We 
further improve upon the method by implement- 
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ing a Monte Carlo loop in which wc performed the 
fit 100 times, each time perturbing the observed 
line flirxes by their measurement uncertainties, as- 
suming a Gaussian distribution. In this way we 
are able to obtain a variance on the DEM at each 
temperature as well as on each elemental abun- 
dance according to the quality of the line measure- 
ments (systematic uncertainties in the calibration 
or the atomic database are, however, still present). 
After a DEM was obtained, wc re-fit the lines us- 
ing the plasma continuum predicted by that DEM. 
This iteration was done three times. We also used 
the synthetic spectrum to adjust the abundance 
and DEM scale factor to match the observed line 
to continuum ratio, and hence the absolute abun- 
dance scale relative to the Solar values of Anders 
& Grevesse (1989). During the iterations, we com- 
pared the model and observed spectra in detail to 
examine various line series and were able to iden- 
tify and include weaker features, and reject some 
features as probable blends or mis-identifications. 

We show the observed and synthetic counts 
spectra in Figures 4a-4d, the DEM in Figure 5, 
and list the abundances in Table 3. The line mea- 
surements and predicted fluxes are listed in Ta- 
ble 2. 

Our method imposes some implicit regulariza- 
tion by smoothing the DEM over a few temper- 
ature grid points. In fitting all elements simul- 
taneously, we can couple the DEM over broader 
temperature ranges and simultaneously derive self 
consistent relative abundances and DEM. This is 
necessary to fit the low temperature regions which 
are not covered by the temperature ranges of the 
emissivities of the iron ions available with HETGS. 
By including fluxes from EUV lines (e.g., Fe ix- 
XVl), we are able to constrain the; low temper- 
ature portion of the DEM (logT < 6.4), which 
N VII and O vii sample in the HETGS spectrum. 
We tested consistency by fltting elements indepen- 
dently or in small groups. While flts were of lower 
quality and limited in temperature range, they 
were consistent, so we performed the final fit over 
all elements, ignoring the density-sensitive He- 
like forbidden and intersystem lines, some strong 
blends (e.g., Ne ix-|-Fe xix 13.4 A), and some ob- 
viously discrepant features which are likely to be 
either mis-identifications, compromised by blends, 
or lines for which emissivities are not accurately 
known. 



Since there arc poorly quantified systematic un- 
certainties in the atomic data (and thus emissivi- 
ties), we repeated the DEM flts with a lower limit 
to the flux uncertainty of 25% (regardless of the 
counts) to globally approximate the atomic data 
systematic errors. The resulting DEM structure 
was the same, but with appropriately larger un- 
certainties. This lends confidence that the DEM 
structure is real and not due to deficiencies in 
atomic data. 

At no point do we formally minimize a binned 
spectrum against a binned synthetic spectrum. 
This had been the only option and the norm for 
low-resolution X-ray spectral modeling with previ- 
ous missions such as ASCA and ROSAT. There are 
many ways for such "global" fltting to fail for high- 
resolution spectra: inaccurate model wavelengths 
can lead to a mis-match between predicted and ob- 
served lines; any spectral region may have features 
missing from the emissivity database; emissivities 
may be inaccurate; or continuum bins can domi- 
nate x^, to list a few pitfalls. For these reasons, 
we prefer a strictly line-based analysis because we 
can more flnely manipulate the features to be fit 
according to their ionic sequences, density sensi- 
tivity, blending, temperature of formation, or any 
other parameter available in the atomic database. 
These techniques have long been in use for UV and 
EUV emission line spectroscopy of coronal plas- 
mas. 

The DEM wc obtain is dominated by two large 
peaks, at logT ~ 6.9 and ~ 7.4. There is a hot 
tail imposed by the presence of Fe xxv, Fe xxvi, 
Ca XX, and Ar XVIII (within large uncertainty, 
since the lines are of low signal to noise ratio) and 
the short wavelength continuum. At logT = 6.2, 
there is a weak peak required by the EUV lines 
from lower ionization states of iron (iX-XVl) which 
provide a much better constraint than does the rel- 
atively low signal-to-noise detection of N Vii (24.8 
A). The overlapping temperature distributions of 
N and the EUV Fe lines serve to better constrain 
the abundance of N. We note that the derived 
DEM is similar to that obtained for the RS CVn 
system HR 1099 by Drake et al. (2001) based on 
Chandra HETG spectra, while the double-peaked 
structure at high temperatures is also reminiscent 
of the emission measure distributions for AR Lac, 
HR 1099, II Peg, and other stars derived from 
EUVE spectra by Griffiths & Jordan (1998) and 
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Sanz-Forcada, Brickhousc & Duproc (2002); the 
cool bump is also seen in some of the emission mea- 
sures derived by the latter authors. It is tempting 
to conclude that the detailed structure in the DEM 
derived in this study and others does indeed reflect 
the true source temperature structure. We cau- 
tion, however, that such apparent structure can 
also arise as a result of errors in the underlying 
atomic data (see, e.g., Kashyap & Drake 1998). 
In particular, errors in ionization equilibria could 
plausibly induce spurious structure since such er- 
rors would be highly correlated with temperature. 
We therefore emphasize that detailed interpreta- 
tion of the DEM structure should proceed with 
caution. 

The very hot portion of the DEM appears to be 
predominantly due to the flares, since the modula- 
tions from Fe XXV and the nearby continuum are 
nearly 90%. Much of the large peak (logT = 7.2- 
7.6) is contributed by flares, having a modulation 
greater than 50%. Below 7.0 there is little mod- 
ulation. This greater variability in the hottest 
plasma component has been seen in other active 
stellar coronae (Bowyer, Drake & Vennes 2000), 
such as II Peg (Huenemoerder, Canizares & Schulz 
2001) during flaring, in Capella (Brickhouse et al. 
2000) in widely separated observations, and in 
the RS CVn stars studied with EUVE by Sanz- 
Forcada, Brickhouse & Dupree (2002). 

4.3. Abundances 

Relative abundances arc determined by the 
DEM fitting procedure, and are strongly corre- 
lated with the DEM solution. An element which 
is isolated in temperature from other ions will be 
degenerate in the DEM normalization and abun- 
dance, since the flux is determined by their prod- 
uct. By fitting all ions simultaneously wc remove 
some of the degeneracy since a series of interme- 
diary strongly overlapping emissivities can couple 
ions which only overlap slightly. By performing a 
Monte Carlo iteration we derive an estimate of the 
range of solutions allowed by measurement uncer- 
tainties. 

The abundances we derived are listed in Ta- 
ble 3, and we plot them against the first ioniza- 
tion potential (FIP) in Figure 6. The abundances 
range from about half to about double the So- 
lar photospheric values, with no simple trend with 
FIP. Neon is about 1.7 times the accepted Solar 



value and three times the iron abundance. This 
neon to iron ratio is a fairly reliable determina- 
tion since the lines form in overlapping tempera- 
ture ranges and so are relatively independent of 
the DEM (assuming they form from ions within 
the same volume). If we average Al and Ca at the 
lowest FIP, wc find them near Solar, and signifi- 
cantly different from the average of Fe, Mg, and Si. 
Al and Ca form in the same temperature range as 
Si, so this difference is also relatively independent 
of the details of the DEM. 

The photospheric abundances of AR Lac 
(Gehren, Ottmann & Reetz 1999) of iron, silicon, 
and magnesium (all low FIP elements) are system- 
atically higher than the coronal values. Taking the 
error- weighted means, we find the relative abun- 
dance to be 0.6 ±0.06 in the coronae, as compared 
to 1.0 ± 0.13 in the photosphere. The average of 
the lowest FIP elements, Al and Ca, is 1.2 ± 0.2, 
which is comparable to the photospheric value. 

Kaastra et al. (1996) and Singh, White & Drake 
(1996) analyzed the same ROSAT and ASCA data 
with different methods, but obtained statistically 
identical results for everything but neon. Their 
values were about half of ours: 0.3 ± 0.02 (again 
averaging over Fe, Si, Ca, and Mg). The latter 
authors gave a good synopsis of the systematic 
uncertainties in low resolution spectral fits and do 
show how changes of 50% can arise from fitting dif- 
ferent spectral regions. Both studies also obtained 
values systematically lower than ours for all other 
elements fit (S, O, N, Ar, Ne). Global fits to 
low-resolution Beppo-SAX spectra of AR Lac by 
Rodono et al. (1999) obtained an average metal 
abundance of 0.66, similar to our values for Fe, Si, 
S, Mg, and O. 

In comparison to previous measurements from 
low rcsohition data, wc arc most discrepant with 
Singh, White & Drake (1996) in the sulfur abun- 
dance, which they claimed was robust even for low 
resolution data, and also differ greatly in Al and 
Ca. We provide a more robust nitrogen abun- 
dance, for which they suspected large calibration 
systematics. The high FIP element abundances 
are near (N, Ar) to above (Ne) Solar photospheric 
values; while there are no photospheric measure- 
ments for these elements for AR Lac, it seems 
likely that they follow the approximately solar val- 
ues for the elements studied by Gehren, Ottmann 
& Reetz (1999). 
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4.4. Line Shapes 

Both stellar components of AR Lac have been 
seen to be X-ray active (Walter, Gibson & Basri 
1983; Siarkowski et al. 1996). We have not been 
able to detect eclipses with our limited phase cov- 
erage, but the spectral resolution and phase cov- 
erage do permit us to search for line profile vari- 
ations. At quadrature, the radial velocity separa- 
tion of the binary components is 230 kms~^ (see 
Gehren, Ottmann & Reetz 1999, for a collation of 
system parameters). The instrumental full- width, 
half-maximum (FWHM) of 0.02 A for MEG and 
0.01 A for HEG yields 316 kms"^ at O viii (19 
A) and 200 kms^i for HEG at Fc xvii (15 A). 
We fit the summed quadrature spectra (OID 6 and 
9) with one or two instrumental profiles (no ther- 
mal or turbulent components included), and with 
a single broadened Gaussian convolved with the 
instrumental response. The lines were not well 
fit by single instrumental profiles. Equally good 
fits could be obtained with either a single Gaus- 
sian with a dispersion of about 0.005-0.01 A, or 
by two Gaussians. The O viii two-Gaussian fit is 
close to the separation expected if both stars are 
active and the activity is localized near their re- 
spective photospheres; the 90% velocity limits are 
at about the instrumental resolution and overlap 
the orbital velocities, -162(-196 : -103) kms"^ 
for the G star, and 123(55 : 197) kms'^ for the K 
star (numbers in parenthesis are 90% confidence 
limits). These values are in reasonable agreement 
with the orbital velocities of —115 and 115 kms"-'^, 
respectively, with a 30 km range within the ob- 
servation. 

To test for broadening in a way which does not 
depend on the calibration of the instrumental pro- 
file, we compared line profiles between the quadra- 
ture and conjunction phases. Both the MEG 
O VIII and HEG Ne-Fe 12A blend were broader 
than the instrumental profile at quadrature, and 
were consistent with the instrumental profile at 
eclipses. We show the profiles and differences in 
Figure 7. 

4.5. Density 

The helium-like triplet lines are well known 
density diagnostics (Gabriel & Jordan 1969, 1973; 

Pradhan & ShuU 1981; Porquct & Dubau 2000). 
The critical densities increase with atomic num- 



ber and arc sensitive over the ranges in log iVg 
of about 10-12 for O vii (AA21.6, 21.8, 22.1, 
logT„a^ = 6.3), 11-13 for Ne ix (AA13.45, 13.55, 
13.70, logT^ax ^ 6.6), and 12 14 for Mg XI 
(AA9.17, 9.23, 9.31, log T^a^: = 6.8), which span 
ranges of interest for coronal plasmas (wavelengths 
refer to the resonance, intcrsystcm, and forbidden 
lines, respectively). The density determination de- 
pends upon flux ratios including the weak intersys- 
tcm lines. Positive detection requires high signal, 
accurate continuum, and resolution of blends. If 
we examine the spectra (Figure 4), we see that the 
continuum is fairly well modeled, so this is not the 
limiting factor in our measurements. The model 
significantly underestimates the Ne ix (Figure 4b) 
and marginally overestimates O vii (Figure 4d) 
forbidden lines. The Mg xi resonance and forbid- 
den lines match well (Figure 4b). 

None of our ratios give a good density con- 
straint: O VII is weak, Ne ix is seriously blended 
with Fe and is listed here only for completeness, 
and the Mg xi intcrsystcm line is weak and pos- 
sibly blended with high-n Ne x H-like transitions. 
The formal ratios indicate logarithmic densities 
(cm~^) of about 10.8 from oxygen, but with 90% 
uncertainties which span the range from 9-12; for 
neon, 11 with 90% uncertainties up to 11.5; for 
magnesium, 12.2 with a 90% upper limit of about 
12.8. Neon and magnesium can be considered to 
define upper limits, but the lower limits are un- 
constrained. The neon intersystem line could be 
up to about 30% iron blends (Ness 2002, Ness et 
al., 2003, in preparation), which would lower the 
upper limit somewhat. 

Hence, wc tentatively conclude that wc have 
detected densities on the order of logA^e ^ H, 
high enough to affect O VII and Ne IX ratios, but 
not Mg XI. There is no reason, however, for the 
values to be identical for different ions, since they 
form at different temperatures. Without better 
measurements of the intcrsystcm lines, we do not 
have data to constrain Nf,{T) more rigorously. 

5. Discussion 

The high resolution X-ray spectrum of AR Lac 
provides a wealth of new information about the 
system. While we have derived more detailed 
emission measure models and abundances than 
possible from low resolution data, there are still 
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inadequacies in the models and significant prob- 
lems to be solved. We must remember that 
there are two stars' coronae involved, and that 
they arc highly variable, making comparison with 
other epochs diSicult. The composite nature is 
not a problem for derivation of emission measures 
since that quantity docs not assume any geometric 
structure. Interpretation of the emission measure, 
however, does require information or assumptions 
about the geometrical distribution of plasma. 

The line widths, being barely resolved by 
HETGS, show that both stars are active, as has 
long been known from the first X-ray observations 
(Walter, Gibson & Basri 1983) as well as from 
spectroscopy in the ultraviolet (Neff et al. 1989; 
Pagano et al. 2001) and optical (Frasca et al. 
2000). Spatial structure can be derived from line 
velocities or occultations. The only conclusive 
modulation is from flares, which have high enough 
frequency and duration (also noted by Pagano 
et al. 2001, for Mg ii) to make it very difficult 
to infer spatial structure. Prior studies suggested 
the G-star to have a compact and uniform chromo- 
sphere, and compact but highly structured corona 
(Pagano et al. 2001, eg). The current lack of mod- 
ulation during primary eclipse egress coupled with 
line broadening at quadrature is consistent with 
this view. The modulation seen by Siarkowski 
et al. (1996) seems to be ruled out; their ASCA 
light curves showed a deep and long primary min- 
imum, which they modeled as extended emitting 
structures between the two stars. Whether such 
extended material is indeed present will require 
better phase coverage at different epochs as well 
as comparison to coronal models in stars which 
should not have binary interaction. The high fre- 
quency and amplitude of flare variability, however, 
may mean that image reconstruction of the corona 
is not feasible (Pease et al. 2002, Drake et al., 2003 
(in preparation)). 

The DEM and abundance model does predict 
the spectrum reasonably well (see Figure 4), but 
there arc still some very large discrepancies for in- 
dividual lines. Some iron line intensities are well 
modeled, but others exhibit significant discrepan- 
cies. For example: Fe xviii A14.208 is significantly 
stronger in the model, while the neighboring blend 
of Fe XVIII A14.256, Fe xx A14.267 is weU modeled 
in the predicted two-to-one ratio. The Fc xvil 
lines, as a series, are also poorly fit: A17 pair 



model is slightly weak, A 16. 78 model is slightly 
strong, while the A15.01 model is nearly twice as 
strong as the data. This latter discrepancy was 
also noted for II Peg (Hucncmoerder, Canizarcs & 
Schulz 2001). Instead of interpreting the apparent 
weakness of the A 15. 01 line as opacity, we suspect 
that the cmissivitics may need revision. The new 
Fe L-shell calculations of Gu (2002) and by Doron 
& Behar (2002) indicate that some ratios may dif- 
fer from earlier calculations by about a factor of 
1.5. We are working to include updated emissivi- 
ties into the synthetic spectrum. 

There could also be deficiencies in the DEM and 
abundance model. We could increase the model 
strength of Fe XVII A15.01 by increasing the DEM 
on the low-temperature side of its emissivity dis- 
tribution (logT ~ 6.5-6.6), so as to minimize the 
affect on Fe xviii. This, however, would have a 
ripple of side affects: the abundance of Mg and 
Ne would have to be reduced, but enhancing the 
cool DEM would also change the ratios of H-like 
Mg XII and Ne x to their He-likc states, since the 
H-like emissivities span the hot peak of the DEM. 
The emissivities of neon, in turn, overlap signifi- 
cantly with those of oxygen, and so forth. We have 
already implicitly optimized the balance of abun- 
dances in the minimization. However, there may 
be local minima, or our solution may be skewed by 
erroneous iron emissivities, or may be too smooth 
to find very sharp temperature structure. Line 
blending is also a significant problem, and will re- 
quire a non-local approach in which isolated lines 
of a blending species are used to predict the con- 
tribution of that species to blended features. We 
have also assumed that each stellar binary com- 
ponent has the same composition, which might be 
reasonable but is unverified. Future efforts will 
address these issues. 

The normalized line flux residuals (x = (data — 
model) /a) are plotted in Figure 8. The scatter is 
quite large, with x^/^^ ~ 3. The Fe xvii lines form 
near logT = 6.7 6.8, and we can see that we have 
lines both stronger and weaker than the model, as 
well as a number which are well fit. In addition to 
the significance of deviations, the bottom panel of 
the figure shows the flux ratio to the model, which 
indicates the percentage deviations. We note that 
the reduced of the spectral energy distribution 
itself is about 1.0; this is misleading, because it 
is dominated by a large number of well modeled 
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continuum bins with relatively low signal-to-noise 
ratios. This illustrates the problem that minimiza- 
tion of a model against a binned spectrum instead 
of against just the extracted line fluxes places less 
statistical weight on spectral lines thought to be 
useful diagnostics, and so devalues the line flux 
information. 

There is at least one other potential source of 
error in the DEM modeling. We have assumed col- 
lisional ionization equilibrium (CIE) , in which col- 
lisional excitation and ionization from the ground 
state, followed by radiative decay, recombination, 
and di-electronic recombination are the dominant 
processes. We have previously argued (Huenemo- 
erder, Canizares & Schulz 2001) that this is ac- 
ceptable even if there is a large flare contribution 
to the flux, since most lines are expected to ther- 
malize quickly (Golub, Hartquist & Quillen 1989; 
Mewe et al. 1985; Doschek et al. 1980). That CIE 
is valid in all regions of the plasma, though, is 
an assumption. If flaring is very frequent, or if 
the apparently quiescent emission were to arise as 
a result of a superposition of many unresolved, 
weaker flares, then the plasma might be driven 
out of CIE. In this regard, improving the model is 
not a trivial proposition, requiring a detailed at- 
tention to each ion series in conjunction with qual- 
ity assessment of model emissivities and more rig- 
orous incorporation of atomic data uncertainties 
(including uncertainties in the ion populations) to 
determine whether or not the data are indeed well- 
described by the standard optically-thin plasma in 
thermal equilibrium. We will accept the current 
DEM and model parameters as the best available, 
given these caveats. 

The abundance trends found in other systems 
show that low PIP species arc sub-Solar, while 
high PIP elements (neon) which have enhanced 
abundance (Huenemoerder, Canizares & Schulz 
2001; Drake et al. 2001; Audard, Giidel & Mewe 
2001; Brinkman et al. 2001; Giidel et al. 2001a,b). 
However, we find that the lowest PIP elements 
have intermediate abundances, as does the high 
PIP element Ar, relative to the other elements. 
AR Lac seems to be more moderate in terms of 
abundance anomalies compared to other active 
stars like II Peg or HR 1099, being less deficient in 
iron, and having a lower neon abundance ratio to 
iron (see Kastner et al. 2002; Huenemoerder 2002, 
for collections of HETGS spectra qualitatively or- 



dered by iron to neon abundance ratio). 

The trend of increasing abundance with PIP is 
opposite what has been seen in the Solar corona 
(see, e.g., Peldman & Laming 2000; Peldman 1992; 
Peldman, Widing & Lund 1990, for reviews). PIP- 
based fractionation is believed to occur in the stel- 
lar chromosphere where low PIP species are pre- 
dominantly ionized while high PIP species remain 
neutral. However, there is as yet no widely ac- 
cepted quantitative explanation of the solar PIP 
effect. The new generation of stellar observations 
from Chandra and XMM-Newton have provided 
new challenges for future models. The abundances 
derived here for AR Lac further complicate the 
picture with enhanced abundances for elements 
with both low and high PIP, relative to interme- 
diate PIP elements. 

The integrated emission measure we obtain for 
the range, logT = 6.0-8.5, of 1.2 x 10^^ cm"^, 
is similar to previous determinations. About one 
third of the emission measure is in the peak at 
logT = 6.9, about half in the peak at 7.4, 15% in 
the hot tail, and about 2% in the cool bump at 
6.2 (note that Pigure 5 plots the DEM integrated 
over intervals of A logT = 0.05). Griflaths & Jor- 
dan (1998) derived emission measures from UV, 
EUV, and X-ray data, and Rodono et al. (1999), 
who analyzed SAX data, found temperature com- 
ponents similar to our peak DEM temperatures, 
and integrated emission measures of about half 
our value. Singh, White & Drake (1996) also de- 
rived similar temperature components and emis- 
sion measures from ROSAT and ASCA spectra, 
with an integrated emission value more similar to 
ours. These diverse observations and analyses over 
several epochs are remarkably similar, and indi- 
cate that the mean activity level of the AR Lac 
coronae is relatively stable. 

If we interpret the emission with a very simple 
geometric model of a semi-toroidal loop of con- 
stant cross section, we can derive a loop height 

as 

(4) 

in which iVioo is the number of identical loops di- 
vided by 100, ao.i is the ratio of loop radius to 
length divided by 0.1, VEM is the volume emis- 
sion measure, N^. is the electron density, and i?* 
is the stellar radius (see Huenemoerder, Canizares 
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& Schulz 2001, for a derivation). If wc attribute 
half the emission to each star and use a density 
of 1 X lO^^cm"^, then the loop heights are 0.08 
and 0.04 (in stellar radii, relative to the G- and 
K-star components, respectively). In other words, 
the loops are compact; even if there are only 10 
loops instead of 100, they only double in height, 
but if density were also ten times lower, then the 
extent becomes significant. If we attribute the hot 
part of the DEM to a single flare loop (since "hot" 
flux is highly modulated), then that loop height is 
about 0.2-0.9 times the radius of the G-star, for 
density values of 0.1 1 x 10^^; this is significant 
compared to the size of the binary system (2-15% 
of the semi-major axis). The simple assumptions 
of semi-toroidal loop intersecting a planar atmo- 
sphere are not valid if loops are large relative to 
the stellar radius, but this heuristic argument sup- 
ports the existence of extended coronal structures 
as has been suggested by several authors (Wal- 
ter, Gibson & Basri 1983; Siarkowski et al. 1996; 
Frasca ct al. 2000; Pagano et al. 2001; Trigilio et al. 
2001). From radio interferometry, Trigilio et al. 
(2001) derived a scale for the "core" component of 
the radio corona of about 0.4 G-star radii, simi- 
lar to our low-density, few-loop case. Their "halo" 
component was somewhat more extended, to sev- 
eral stellar radii. It is, however, difficult to make 
a meaningful comparison between the X-ray and 
radio extents due to the very different model as- 
sumptions and the very uncertain X-ray plasma 
density determination, and because the radio and 
X-ray emission to not necessarily originate from 
the same plasma. 

In contrast to the spectral behavior, the light 
curves at various epochs can be quite different. 
The SAX light curve (Rodono et al. 1999) showed 
much structure, with a rotational modulation, a 
short eclipse, and frequent and short flares. In 
contrast, the ASCA light curve (Siarkowski et al. 
1996) was relatively smooth and showed broad and 
deep primary eclipse and broad, shallow secondary 
eclipse. Rodono et al. (1999) smiimarize and com- 
pare many of the X-ray observations. Some as- 
pects of the data are quite confusing and challeng- 
ing to explain. For example, the SAX data show a 
narrow dip during primary eclipse, beginning after 
second contact (G-star photosphere fully eclipsed), 
and returning to the higher level before third con- 
tact (G-star photosphere fully egressed). This is 



not consistent with their interpretation that the 
G-star corona is compact, since the feature would 
have to be trailing the G-star to make ingress 
constant in flux, yet leading on egress. Instead, 
there must be other emerging or erupting struc- 
tures which could be on either star. Our highly 
non-repeatable light curves show the extreme diffi- 
culty of interpreting variability as rotational mod- 
ulation or eclipses of stable structures, as was also 
noted by Pease et al. (2002) for EUV variability. 

The flux of AR Lac integrated from 1.7-30 A 
is 3.7 X 10~^^ ergscm~^s~^. For a distance of 42 
pc, the band luminosity is 7.9 x 10^°ergss~^. This 
also represents the time-averaged value over the 
observation, which is strongly biased to the longer 
quadrature integrations at lower count rate. For 
short times, the luminosity can be several times 
higher; the flare of OID 7 emitted about 10^^ ergs. 
This is a typical size for RS CVn binaries (e.g. 
Maggio et al. 2000; Huenemoerder, Canizares & 
Schulz 2001) but there have been flares two orders 
of magnitude larger (Ottmann & Schmitt 1994; 
Kuerster & Schmitt 1996). Perhaps the flare we 
missed in X-rays was of such a magnitude, given 
the large X-ray to EUV enhancement ratio we see 
near orbit 2.4 (see Figure 2, top panel). 

Detailed independent modeling of flare and 
non-flare states, while scientiflcally interesting in 
terms of temperature and abundance structure, 
unfortunately requires better statistics than we 
have in these data. The larger X-ray flare of OID 
7, while significant in the light curve, is relatively 
short for DEM analysis — much of the count- 
rate change is from the continuum. The EUVE 
data, even binned over the entire observation, 
serve primarily to constrain the low temperature 
emission since the signal is so miich weaker than 
the HETGS data where lines overlap in formation 
temperature. We also looked for spectral differ- 
ences near phase 0.15 between OID 6 and 9, but 
signal was not adequate, even for summed line 
fluxes. The line modulation shown in Figure 3 
is our low-signal proxy for the detailed model. 
Time-dependent DEM and abundance modeling 
of flares will have to await larger events. 

We might expect large flares to result in line 
shifts or excess broadening due to velocity fields. 
In Figure 7 we show quadrature and conjunction 
profiles for two lines, O viii and Ne x, which 
are near the highest resolutions obtained with 
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HETGS. These features, however, form at rel- 
atively low temperatures and are relatively un- 
modulated by the flares (see Figure 3). We have 
attributed the broadening in O viii to orbital ef- 
fects. Ne X is interesting because the conjunc- 
tion profile shows a marginal blueshift, and these 
phases were most affected by flares. Ayres et al. 
(2001) tentatively detected a transient blueshift in 
Ne X in HR 1099, coincident with a UV flare; there 
was no corresponding X-ray flare detected. We 
likewise find the effect inconclusive due to conflict- 
ing information and marginal data quality. How- 
ever, the tentative result is intriguing and should 
foster further observations and more sophisticated 
analyses. 

6. Conclusions 

We have presented an analysis of high reso- 
lution X-ray and EUV spectra and photometry 

of the bright eclipsing RS CVn system AR Lac. 
While our results are qualitatively similar to those 
of earlier studies, the high resolution X-ray line 
spectra have allowed us to obtain a significantly 
more detailed glimpse of the coronal temperature 
structure and abundances. 

The X-ray and EUV spectra show that the coro- 
nae of this system are characterized by compli- 
cated structure and variability. The overall activ- 
ity level is similar to that determined from previ- 
ous observations. The hot portion of the DEM is 
strongly modulated by flares, whereas the cooler 
portion appears relatively quiescent. The flare 
modulation is either frequent and large enough to 
hide eclipses, or the dominant coronal structures 
are of polar origin and are not rotationally mod- 
ulated. The steadily emitting structures could be 
compact if the density estimate of 1 x 10^^ cm"'^ is 
accurate, whereas if the emission from the X-ray 
flare is from a single loop at this density, then it 
would be significantly extended. This ambiguity 
in interpretation illustrates the importance of reli- 
able plasma density estimates and underscores the 
difficulty in obtaining such estimates, even from 
high-quality Chandra HETG observations. Given 
the high resolution of the HETGS, we have deter- 
mined that both stellar components of the system 
contribute significantly to the X-ray emission. The 
quadrature line profiles are consistent with the K- 
and G-stars being equally bright. 



Coronal abundances show some similarities 
with those found for other RS CVn stars but differ 
substantially in detail: the highest and lowest FIP 
species both appear to have higher abundances rel- 
ative to intermediate FIP ions. For AR Lac, the 
intermediate FIP abundances are below measured 
photospheric values. 
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Table 1 
Observational Information. 



OID 


Start Date & Time 


Exposure** 


Phase^ 


6 


2000-09-llT22:05:13 


32.7 


0.14-0.34 


7 


2000-09-16Tll:48:03 


7.7 


0.48-0.52 


8 


2000-09-15T14: 18:04 


9.5 


0.98-0.02 


9 


2000-09-17T20:02:12 


32.5 


0.15-0.35 


10 


2000-09-20T09:28:ll 


7.5 


0.42-0.46 


11 


2000-09-19X14:31:14 


7.5 


0.02-0.06 



'^Exposure time is in ks. 

''Phase is computed for the exposure start and stop 
times using the cphcmcris of Ferryman et al. (1997). At 
phase 0.0, the G-star is totally eclipsed. 
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Table 2 
Emission Line Data. 



Mnemonic^ 


Ion 


log 




At = 


Ao 


d 




fi 


■= 




/ 


t 


St 


g 




Fe26HLa 


Fe XXVI 


8 


1 


1 


781 


1.782 (1 


56 


2) 


3.48 


(4 


43) 


4 


40 


-0 


92 


-0.2 


Fe25HeLa 


Fe XXV 


7 


8 


1 


861 


1.859 


( 1 


9) 


23.70 


(7 


14) 


25 


95 


-2 


25 


-0.3 


Ca20HLa 


Ca XX 


7 


8 


3 


021 


3.020 (1 


10 


5) 


1.97 


(1 


19) 


1 


90 





06 


0.1 


CalQHeLaB 


Ca XIX 


7 


5 


3 


198 


3.192 ( 


18 


1) 


8.10 


(2 


74) 


8 


37 


-0 


27 


-0.1 


ArlTHcLb 


Ar XVII 


7 


4 


3 


365 


3.365 


( 4 


3) 


3.63 


(1 


89) 


1 


16 


2 


47 


1.3 


ArlSHLa 


Ar XVIII 


7 


7 


3 


734 


3.732 


( 4 


9) 


4.44 


(2 


08) 


5 


57 


-1 


13 


-0.5 


ArlTHoLar 


Ar XVII 


7 


4 


3 


949 


3.950 


( 2 


4) 


9.32 


(2 


49) 


8 


91 





41 


0.2 


ArlTHoLai 


Ar XVII 


7 


3 


3 


968 


3.969 


( 5 


6) 


5.49 


(2 


36) 


2 


32 


3 


17 


1.3 


S16HLb,Arl7HcLaP 


Ar XVII 


7 


5 


3 


992 


3.994 


( 3 


5) 


6.10 


(2 


19) 


6 


12 


-0 


02 


-0.0 


S16HLa 


S XVI 


7 


6 




730 


4.730 


( 


9) 


23.60 


(3 


42) 


26 


05 


-2 


45 


-0.7 


S15HeLar 


S XV 


7 


2 


5 


039 


5.039 


( 1 


0) 


24.85 


(3 


56) 


26 


65 


-1 


80 


-0.5 


SlSHeLai 


S XV 


7 


2 


5 


065 


5.065 ( 


13 


8) 


6.45 


(2 


80) 


5 


83 





61 


0.2 


S15HeLaf 


S XV 


7 


2 


5 


102 


5.100 


( 2 


4) 


13.59 


(3 


13) 


8 


78 


4 


81 


1.5 


Sil4HLb 


Si XIV 


7 


4 


5 


217 


5.219 


( 2 


2) 


12.83 


(3 


30) 


8 


81 


4 


02 


1.2 


Sil4HLa 


Si XIV 


7 


4 


g 


183 


6.181 


( 


4) 


68.89 


(2 


94) 


62 


34 


6 


55 


2.2 


SilSHeLar 


Si XIII 


7 





Q 


648 


6.648 


( 


5) 


44.03 


(2 


28) 


51 


71 


-7 


68 


-3.4 


Sil3HeLai 


Si XIII 


7 





g 


687 


6.686 


( 1 


5) 


8.49 


(1 


53) 


9 


54 


-1 


05 


-0.7 


Sil3HeLaf 


Si XIII 


7 





g 


740 


6.740 


( 


3) 


30.20 


(1 


99) 


20 


64 


9 


56 


4.8 


Mgl2HLb 


Mg XII 


7 


2 


7 


106 


7.105 


( 


7) 


13.47 


(1 


47) 


14 


14 


-0 


67 


-0.5 


Fe24w7.17 


Fe XXIV 


7 


4 


'J 


169 


7.170 


( 


7) 


13.84 


(1 


45) 


4 


33 


9 


51 


6.6 


AllSHLaj 


Al XIII 


7 


3 


I 


1 TO 


7.170 


( 


5) 


11.03 


(0 


94) 


12 


00 


-0 


97 


-1.0 


A112HeLar 


Al XII 


7 





1 


lot 


7.758 


( 1 


9) 


7.84 


( 1 


35) 


5 


51 


2 


33 


1.7 


MgllHeLb 


Mg XI 


6 


9 


7 


850 


7.849 


( 2 


6) 


6.00 


( 1 


36) 


7 


66 


-1 


67 


-1.2 


A112HcLaf 


Al XII 


6 


9 


7 


872 


7.869 


( 4 


5) 


4.71 


( 1 


35) 


4 


53 





18 


0.1 


Fo24w7.99 


Fc XXIV 


7 


4 


7 


991 


7.989 


( 1 


6) 


9.48 


( 1 


31) 


13 


82 


-4 


34 


-3.3 


Fc24w8.23 


Fe XXIV 


7 


4 


8 


233 


8.283 


( 4 


2) 


3.24 


( 1 


32) 


4 


96 


-1 


71 


-1.3 


Fe24w8.28 


Fe XXIV 


7 


4 


8 


285 


8.304 


( 2 


5) 


6.57 


( 1 


74) 


1 


85 


4 


72 


2.7 


Fe23w8.30 


Fe XXIII 


7 


2 


8 


304 


8.318 


( 6 


5) 


8.12 


( 1 


95) 


8 


36 


-0 


24 


-0.1 


Fe24w8.32 


Fe XXIV 


7 


4 


o 
o 


'^^fi 
o±o 


8.237 


( 1 


8) 


5.41 


( 1 


30) 


10 


09 


-4 


68 


-3.6 


Mgl2HLa 


Mg XII 


7 


2 


g 


422 


8.420 


( 


3) 


98.76 


( 3 


20) 


102 


70 


-3 


94 


-1.2 


Fe23w8.81 


Fe XXIII 


7 


2 


8 


815 


8.814 


( 1 


7) 


8.34 


( 1 


47) 


8 


38 


-0 


04 


-0.0 


Fe22w8.97 


Fe XXII 


7 


1 


8 


975 


8.975 


( 1 


3) 


9.95 


( 1 


64) 


7 


57 


2 


38 


1.4 


MgllHoLar 


Mg XI 


6 


8 


9 


169 


9.168 


( 


5) 


57.25 


( 3 


01) 


59 


28 


-2 


03 


-0.7 


MgllHcLai 


Mg XI 


6 


8 


9 


230 


9.231 


( 2 


1) 


14.32 


( 2 


06) 


8 


99 


5 


33 


2.6 


MgllHcLaf 


Mg XI 


6 


8 


9 


314 


9.313 


( 


7) 


28.47 


( 2 


07) 


26 


81 


1 


66 


0.8 


NclOHLg 


No X 


7 





9 


708 


9.710 


( 1 


2) 


45.64 


( 6 


46) 


28 


24 


17 


40 


2.7 


NelOHLb 


Ne X 


7 





10 


239 


10.238 


( 


4) 


85.51 


( 3 


80) 


89 


03 


-3 


52 


-0.9 


Fe24wl0.62 


Fe XXIV 


7 


4 


10 


619 


10.620 


( 


9) 


61.70 


( 4 


27) 


65 


97 


-4 


27 


-1.0 


Fe24wl0.66 


Fe XXIV 


7 


4 


10 


663 


10.661 


' 1 


9) 


34.56 


( 3 


54) 


34 


61 


-0 


05 


-0.0 


Fel7wl0.77 


Fe XVII 


6 


8 


10 


770 


10.768 


( 2 


5) 


10.41 


[ 2 


09) 


9 


08 


1 


32 


0.6 


Fel9wl0.82 


Fe XIX 


6 


9 


10 


816 


10.814 


( 2 


5) 


11.44 


' 2 


16) 


11 


96 


-0 


52 


-0.2 


Fe23wl0.98 


Fe XXIII 


7 


2 


10 


981 


10.983 


( 1 


2) 


44.19 


( 7 


54) 


44 


10 





09 


0.0 


Fc23wll.02 


Fe XXIII 


7 


2 


11 


019 


11.015 


( 3 


8) 


24.62 


( 7 


90) 


28 


89 


-4 


27 


-0.5 


Fo24k'11.03 


Fe XXIV 


7 


4 


11 


029 


11.031 


( 7 


8) 


54.16 


( 9 


48) 


42 


19 


11 


97 


1.3 


Fo24wll.l8 


Fe XXIV 


7 


4 


11 


176 


11.173 


( 


6) 


72.06 


( 4 


36) 


76 


20 


-4 


14 


-1.0 


r eiow±± .Jo 


re A. V ill 


6 


8 


11 


326 


1 1.330 





9) 


18.26 


2 


D4 ) 


18 


61 


-0 


35 


-0.1 


Fel8wll.53 


Fe XVIII 


6 


8 


11 


527 


11.529 


( 2 


2) 


19.78 


( 3 


11) 


12 


78 


7 


00 


2.3 


Ne9HeLb 


Ne IX 


6 


6 


11 


544 


11.549 


' 1 


2) 


28.16 


( 3 


37) 


20 


24 


7 


92 


2.3 


Fe23wll.74 


Fe XXIII 


7 


2 


11 


736 


11.738 


( 1 


1) 


91.55 


( 4 


78) 


91 


19 





36 


0.1 


Fe22wll.77 


Fe XXII 


7 


1 


11 


770 


11.771 


( 1 


1) 


74.99 


( 4 


45) 


70 


06 


4 


93 


1.1 


Fel7wl2.12'' 


Fe XVII 


6 


7 


12 


124 


12.130 


( 


9) 


251.30 ( 


33 


48) 


50 


17 


201 


13 


6.0 


NelOHLa'' 


Ne X 


6 


9 


12 


135 


12.139 


( 3 


8) 


383.89 ( 


33 


00) 


633 


70 


-249 


81 


-7.6 


Fc23wl2.16 


Fe XXIII 


7 


2 


12 


161 


12.155 


( 1 


9) 


56.61 ( 


10 


70) 


50 


26 


6 


35 


0.6 


Fol7wl2.27 


Fe XVII 


6 


7 


12 


266 


12.261 


( 3 


4) 


27.94 


( 9 


27) 


45 


09 


-17 


16 


-1.9 


Fc21wl2.28 


Fe XXI 


7 


1 


12 


284 


12.285 


( 


6) 


137.49 ( 


15 


15) 


135 


90 


1 


59 


0.1 


Fc20wl2.58 


Fe XX 


7 





12 


576 


12.574 


( 


0) 


21.54 


( 3 


34) 


22 


67 


-1 


13 


-0.3 


Fe22wl2.75 


Fe XXII 


7 


1 


12 


754 


12.751 


( 1 


9) 


23.90 


( 3 


63) 


25 


34 


-1 


44 


-0.4 


Ne9HeLar 


Ne IX 


6 


6 


13 


447 


13.445 


( 


6) 


132.67 ( 


11 


30) 


144 


90 


-12 


23 


-1.1 


Fel9wl3.50 


Fe XIX 


6 


9 


13 


497 


13.504 


( 


9) 


69.49 


( 7 


51) 


44 


98 


24 


50 


3.3 


Fo19k'13.52 


Fc XIX 


6 


9 


13 


518 


13.520 


( 


6) 


81.26 


( 7 


61) 


99 


15 


-17 


89 


-2.4 


No9HoLai 


Nc IX 


6 


6 


13 


552 


13.554 


( 1 


9) 


36.59 


(4 


81) 


19 


55 


17 


04 


3.5 


Nc9HcLaf 


Nc IX 


6 


6 


13 


699 


13.699 


( 


6) 


110.30 


( 7 


00) 


65 


60 


44 


70 


6.4 


Fol8wl4.26 


Fe XVIII 


6 


8 


14 


256 


14.205 


( 


3) 


119.10 


( 8 


24) 


40 


82 


78 


28 


9.5 


Fe20wl4.27 


Fe XX 


7 





14 


267 


14.259 


( 1 


6) 


48.42 


(6 


15) 


26 


94 


21 


48 


3.5 


Fel8wl4.53 


Fe XVIII 


6 


8 


14 


534 


14.540 


( 1 


3) 


52.72 


(5 


78) 


41 


11 


11 


61 


2.0 
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Table 2 — Continued 



■ f 

iVl 11 c mo 111 c 




Ion 


log J 


)^ 


t 






Jl 




Jt 


t 








h 


OSHLd 


O VIII 


6.7 


14 


821 


14.820 ( 


0) 


21.59 ( 4 


95) 


17 


90 


3 


69 





7 


Fcl7wl5.01 


Fo XVII 


6.7 


15 


014 


15.012 ( 


6) 


271.60 ( 12 


64) 


441 


60 


-170 


00 


-13 


4 


Fcl9wl5.08 


Fo XIX 


6.9 


15 


079 


15.079 ( 1 


3) 


44.58 ( 6 


18) 


33 


30 


11 


28 


1 


8 


08HLg' 


O VIII 


6.7 


15 


176 


15.180 ( 2 


2) 


82.16 ( 8 


92) 


40 


88 


41 


28 


4 


6 


Fel7wl5.26 


Fe XVII 


6.7 


15 


261 


15.261 ( 1 


2) 


116.33 ( 9 


19) 


124 


40 


-8 


07 


-0 


9 


Fel8wl5.62 


Fe XVIII 


6.8 


15 


625 


15.625 ( 1 


9) 


42.85 ( 6 


16) 


55 


81 


-12 


96 


-2 


1 


Fel8wl5.82 


Fe XVIII 


6.8 


15 


824 


15.826 ( 2 


2) 


32.92 ( 5 


94) 


33 


99 


-1 


06 


-0 


2 


Fol8wl5.87 


Fc XVIII 


6.8 


15 


870 


15.869 ( 1 


9) 


33.99 ( 6 


10) 


18 


03 


15 


96 


2 


6 


08HLb"' 


O VIII 


6.7 


16 


006 


16.005 ( 


9) 


130.71 ( 26 


12) 


127 


10 


3 


61 





1 


Fcl8wl6.07 


Fc XVIII 


6.8 


16 


071 


16.072 ( 


9) 


99.86 ( 9 


51) 


77 


00 


22 


86 


2 


4 


Fcl9wl6.11 


Fc XIX 


6.9 


16 


110 


16.106 ( 3 


4) 


25.85 ( 6 


84) 


43 


53 


-17 


68 


-2 


6 


Fcl8wl6.16 


Fc XVIII 


6.8 


16 


159 


16.168 ( 3 


7) 


21.99 ( 6 


20) 


31 


16 


-9 


17 


-1 




Fcl7wl6.78 


Fc XVII 


6.7 


16 


780 


16.773 ( 1 


2) 


168.16 ( 12 


85) 


193 


60 


-25 


44 


-2 





Fel7wl7.05 


Fe XVII 


6.7 


17 


051 


17.049 ( 


6) 


236.37 ( 14 


89) 


232 


30 


4 


07 





3 


Fel7wl7.10 


Fe XVII 


6.7 


17 


096 


17.094 ( 


6) 


222.20 ( 14 


59) 


210 


10 


12 


10 





8 


Fel8wl7.62 


Fe XVIII 


6.8 


17 


623 


17.619 ( 2 


5) 


37.17 ( 8 


06) 


55 


56 


-18 


39 


-2 


3 


07HeLb 


O VII 


6.4 


18 


627 


18.625 ( 7 


5) 


25.34 ( 8 


67) 


15 


36 


9 


98 


1 


2 


08HLa 


O VIII 


6.7 


18 


970 


18.968 ( 2 


0) 


810.96 ( 87 


37) 


858 


80 


-47 


84 


-0 


5 


07HoLar 


O VII 


6.3 


21 


602 


21.601 ( 1 


9) 


125.43 ( 22 


71) 


128 


90 


-3 


47 


-0 


2 


07HoLai 


O VII 


6.3 


21 


802 


21.804 (156 


2) 


23.69 ( 15 


62) 


18 


70 


5 


00 





3 


07HcLaf 


O VII 


6.3 


22 


098 


22.095 ( 11 


3) 


43.55 ( 19 


63) 


76 


97 


-33 


42 


-1 


7 


N7HLa 


N VII 


6.5 


24 


782 


24.780 ( 3 


1) 


132.88 ( 25 


49) 


137 


20 


-4 


32 


-0 


2 


Fcl8w93.92 


Fe 


XVIII 


6.8 


93 


923 


93.920 ( 3 


0) 


661.91 ( 118 


20) 


711 


70 


-49 


79 


-0 


4 


Fel9wl01.55 


Fe 


XIX 


6.9 


101 


550 


101.510 ( 40 


0) 


259.18 ( 78 


54) 


235 


00 


24 


18 





3 


Fel8wl03.94 


Fe 


XVIII 


6.8 


103 


937 


103.940 ( -3 


0) 


250.63 ( 86 


43) 


260 


50 


-9 


87 


-0 


1 


Fcl9wl08.37 


Fe 


XIX 


6.9 


108 


370 


108.230 (140 


0) 


796.49 ( 124 


45) 


900 


40 


-103 


91 


-0 


8 


Fol9wl09.97 


Fc 


XIX 


6.9 


109 


970 


109.790 (180 


0) 


217.01 ( 74 


83) 


120 


20 


96 


81 


1 


3 


Fc20wll0.63 


Fc 


XX 


7.0 


110 


630 


110.830 (-200 


0) 


181.00 ( 78 


69) 


42 


81 


138 


19 


1 


8 


Fc22wll4.41 


Fc 


XXII 


7.1 


114 


410 


114.430 (-20 


0) 


339.81 ( 94 


39) 


358 


90 


-19 


09 


-0 


2 


Fe20wll8.66 


Fe 


XX 


7.0 


118 


660 


118.520 (140 


0) 


284.33 ( 98 


05) 


528 


10 


-243 


77 


-2 


5 


Fel9wl20.00 


Fe 


XIX 


6.9 


120 


000 


119.810 (190 


0) 


177.00 ( 93 


16) 


243 


60 


-66 


60 


-0 


7 


Fe20wl21.83 


Fe 


XX 


7.0 


121 


830 


121.750 ( 80 


0) 


553.19 ( 134 


93) 


1029 


00 


-475 


81 


-3 


5 


Fe21wl28.73 


Fe 


XXI 


7.0 


128 


730 


128.560 (170 


0) 


1445.10 ( 190 


13) 


1520 


00 


-74 


90 


-0 


4 


Fe22wl35.78 


Fe 


XXII 


7.1 


135 


780 


135.590 (190 


0) 


1168.80 ( 201 


53) 


799 


30 


369 


50 


1 


8 


Fc21wl42.16 


Fe 


XXI 


7.0 


142 


215 


142.380 (-165 


0) 


437.97 ( 199 


07) 


139 


70 


298 


27 


1 


5 


Fc9wl71.07 


Fc 


IX 


5.8 


171 


073 


170.950 (123 


0) 


1240.10 ( 539 


18) 


1377 
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Note. — Values given in parentheses are the one standard deviation uncertainties on the preceding quantity. Lines with 
wavelength greater than 30 A are EUVE data. 

^Thc mnemonic is a convenience for uniquely naming each feature. It is comprised of the clement and ion (in Arabic 
numerals) followed by a string indicating a wavelength and the wavelength (e.g., wl6.78), or a code for the hydrogen- like 
("H'') and helium-like "He" series, "L" for Lyman transition, one of "a", "b", "g", "d", or "e" for series lines a, /3, 7, 5, e, 
and "r" , "i" , or "f" for resonance, intersystem, or forbidden lines. 

"^Average logarithmic temperature [Kelvins] of formation, defined as the first moment of the emissivity distribution. 

"^Theoretical wavelengths of identification (from APED), in A. If the line is a multiplet, we give the wavelength of the 
strongest component. 

'^Measured wavelength, in A (uncertainty is in mA). 

^Emitted source line flux is 10~^ times the tabulated value in [phot cm~^ s~^]. 
^Model line flux is 10~® times the tabulated value in [photcm~^ 
®Line flux residual, Sf = fo — ft. 
""Sx^Uo- ft)/^o. 

XVI 3.991 A is blended with Ar xvii 3.994 A in about equal strengths for the model DEM. 

JAl xiii 1.172 is blended with Fe xxiv 7.169. We have adjusted the Al flux to account for a 25% contribution of Fe for 

the assumed DEM. 

^Ne X 12.132 A is blended with Fc xvn 12.124 A; in the HEG spectrum, the lines are barely resolved, and expected to 
be about 10% as strong. Their sum has 5% — —0.3. 

VIII 15.18 A is blended with Fe xix in about equal strengths. 

"^O VIII 16.01 A is l>lcndcd with Fe xvill; Fe XVIII 15.62 A has a similar theoretical emissivity distribution to Fe XVIII 
16.01, so we have subtracted the latter's flux from the me^yred O viii to approximate the net O viii flux. 



Table 3 
Elemental Abundances 



Element 




Corona*^ 


Photosphere'^ 


N 


14.53 


1.2 (0.3) 







13.62 


0.6 (0.1) 




Ne 


21.56 


1.6 (0.3) 




Mg 


7.65 


0.7 (0.1) 


1.0 (0.2) 


Al 


5.99 


1.3 (0.3) 




Si 


8.15 


0.6 (0.1) 


1.3 (0.9) 


S 


10.36 


0.6 (0.1) 




Ar 


15.75 


1.0 (0.2) 




Ca 


6.11 


1.0 (0.3) 


1.3 (0.4) 


Fe 


7.87 


0.5 (0.1) 


1.0 (0.2) 



^FIP is the first ionization potential in eV. 

Coronal abundances are fractional relative to the 
Solar abundances of Anders & Grevesse (1989). Num- 
bers in parentheses are one standard deviation uncer- 
tainties determined from Monte Carlo fitting. 

'^Photospheric abundances from Gehren, Ottmann 
& Reetz (1999). 
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Fig. 1 . — This is the energy flux density spectrum 
for the summed HEG and MEG first orders for 
all observations, totalling 97 ks. The spectrum 
has been smoothed by convolving it with a Gaus- 
sian kernel of dispersion 0.005A. Significant lines 
have been marked. The strength of the short- 
wavelength continuum (< 5 A) and the number 
of high-excitation iron lines are indicative of the 
presence of very hot plasma up to about 100 MK, 
while presence of N vii, O vii, O VIII, and Ne IX 
indicates cooler plasma, near 3 MK. 

Fig. 2. — We show the light curve over the six ob- 
servations in several forms. The top panel shows 
the count rate against a binary system orbit num- 
ber, with zero defined as a primary conjunction 
(HJD 2448501.1232). Given the orbital period 
of 1.98318 days, an orbit is 171 ks. Both MEG 
and HEG orders —3 to -|-3 (excluding the zero 
order, which is saturated) were summed into 0.5 
ks bins over the wavelength range of 1.7-25 A. 
Data are labeled by their OID. The EUVE hght 
curve is plotted with small open squares (colored 
green in the electronic edition), and runs from 
about orbit 1.6 to 3.0, and is in 0.6 ks bins. The 
EUVE data's error bars, which have been omit- 
ted for clarity, are all less than 15%. Note the 
small EUV flux enhancement during the flare in 
OID 7. The next panel shows the same data, 
but phase-folded, and only over the phase inter- 
val observed by Chandra. The dotted curve (solid 
orange in the electronic edition) is a simple oc- 
cultation model for uniform disks of equal surface 
brightness to show where photospheric eclipses oc- 
cur, arbitrarily scaled to exaggerate the eclipses. 
The primary eclipse (phase = 0.0) is total, and 
the secondary is annular. The EUVE data are 
the sparse set of squares running along the bot- 
tom (colored green in the electronic edition). The 
bottom pair of graphs show the light curves for 
narrow bands around Si XIV (6.16-6.22 A), which 
is strongly modulated (upper pane), and Si XIII 
(6.62-6.78 A), which is only weakly modulated 
(bottom pane). Each has had a nearby continuum 
band rate scaled and subtracted. Different orbits 
at the same phase are distinguished by the error- 
bar line style (solid blue in the electronic edition). 

Fig. 3. — We computed the modulation of fea- 
tures from their count-rate, defined as the {max — 
min)/{max + min). Lines are labeled by element 



and ion (in Arabic notation). These arc plotted 
against their temperature of maximum emissivity. 
The hotter plasma is more strongly modulated, 
which is consistent with variations being due to 
flares, rather than a change in volume which would 
affect all lines equally. 

Fig. 4.— These arc the MEG and HEG counts 
spectra and the folded DEM and abundance model 
for the summed -|-1 and —1 orders. Both data and 
model have been smoothed by a Gaussian convolu- 
tion with dispersions of 0.005 A (MEG) and 0.0025 
A (HEG). The model is plotted with a dashed line 
(solid red in the electronic edition). The MEG 
spectrum is usually the upper trace, except near 
Fe XXV, 1.85 A, where the HEG has more effec- 
tive area. HEG has about twice the resolution of 
MEG. Some model lines expected in each region 
have been labeled. Some regions of the spectrum 
are well matched by the model. In others, large 
mismatches arc obvious. The spectra are shown in 
four parts, a-d, each of which has four panels cov- 
ering about 1.5 A. Part a spans 1.7-7.7 A, part b 
7.7-13.7 A, part c 13.7-19.7, and part d 19.7-25.7 
A. 

Fig. 4b.— See 4 
Fig. 4c. — See 4 
Fig. 4d.— See 4 

Fig. 5. — We show the differential emission mea- 
sure integrated in bins of size AlogT = 0.05 for 
the DEM fit to combined observations. The solid 
line is the mean of the 100 Monte-Carlo fits, and 
the dashed lines arc the one standard deviation 
boundaries. The tail above about logT = 7.6 
was manually adjusted (within the one sigma lim- 
its) so that the synthetic spectrum's short wave- 
length continuum better matched the data, since 
this wasn't well constrained by the fit. The peak 
at 6.2 is constrained by contemporaneous EUVE 
data. The range including the hot peak at 7.3 and 
hotter tail is strongly modulated by flares. The 
integral over the plotted temperature range yields 
a volume emission measure of 1.2 x 10^'*cm~^. 

Fig. 6. — The abundances relative to Solar are 
plotted against their first ionization potential. Re- 
sults from the current work are marked with 
squares, and the uncertainties are one standard 



19 



deviation from the Monte Carlo fit. Results from 
low-resolution X-ray determinations are also plot- 
ted (diamonds and asterisks, with dotted error- 
bars) as well as the AR Lac photosphcric values 
(triangles) and scaled Solar coronal abundances 
(dashed line); the symbol legend indicates the first 
author. 

Fig. 7. — In this figure we compare the normalized 
flux profiles at quadrature and conjunction phases. 
Normalization is done within the wavelength band 
by subtracting the minimum and then dividing by 
the maximum. The solid line profile is from the 
conjunction phases when the stars' line-of-sight ve- 
locity difference is zero, the dashed quadrature, 
and the lowest dash-dot line is the difference. The 
left panel is from MEG and the right is HEG (at 
the same wavelength, HEG has twice the resolu- 
tion of the MEG). Both broadening and shifts are 
apparent. 

Fig. 8. — In this figure we plot the line fiux residu- 
als normalized by their measurement uncertainty 
(upper graph) and the ratio of measured line flux 
to model fiux (bottom graph). The abscissa is 
the first moment of the line's emissivity distribu- 
tion, rather than the temperature of maximum 
emissivity. This is a conservative bias for the 
hydrogen-like ions, since they have a long tail on 
the high-tempcraturc side they move to slightly 
higher values than their peak. In the both plots, 
each ion is denoted by the first character in its 
element abbreviation with some use of lower-case 
for duplicate characters (see the legend). While 
many features lie within ztlcr of zero, many lie 
outside two and 3a, in excess of what would be 
expected statistically for Gaussian measurement 
errors alone. Hence, there is a deficiency in the 
DEM and abundance model, errors in the atomic 
data, mis-identified and blended lines, or likely, 
some combination of these. 



20 



> 

X 
X 



>«< 



xx: 



X 

c 



?C 5*^ ^■•"k^ 

>»=X Xx>«< X>S< X Si X 

^<L) 0) <u CUD d) 0^ 0) OJ 

ZI^ULi- b_ZL Li_ 'HZ: Li_ 



» 
XX 



» > 

XX X 



» 
oo 



ODO 



10 



15 

Wovelenath FAl 



20 



25 



Fig. 1. — This is the energy flux density spectrum for the summed HEG and MEG first orders for all 
observations, totalling 97 ks. The spectrum has been smoothed by convolving it with a Gaussian kernel 
of dispersion 0.005A. Significant lines have been marked. The strength of the short-wavelength continuum 
(< 5 A) and the number of high-excitation iron lines are indicative of the presence of very hot plasma up to 
about 100 MK, while presence of N vii, O vii, O viii, and Ne ix indicates cooler plasma, near 3 MK. 
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Fig. 3. — We computed the modulation of features from their count-rate, defined as the {m,ax — min) / (max + 
min). Lines are labeled by element and ion (in Arabic notation). These are plotted against their temperature 
of maximum emissivity. The hotter plasma is more strongly modulated, which is consistent with variations 
being due to flares, rather than a change in volume which would affect all lines equally. 
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Fig. 5. — Wc show the differential emission measure integrated in bins of size AlogT = 0.05 for the DEM 
fit to combined observations. The sohd hne is the mean of the 100 Monte-Carlo fits, and the dashed lines 
are the one standard deviation boundaries. The tail above about logT = 7.6 was manually adjusted (within 
the one sigma limits) so that the synthetic spectrum's short wavelength continuum better matched the data, 
since this wasn't well constrained by the fit. The peak at 6.2 is constrained by contemporaneous EUVE 
data. The range including the hot peak at 7.3 and hotter tail is strongly modulated by flares. The integral 
over the plotted temperature range yields a volume emission measure of 1.2 x lO^^cm"^. 
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Fig. 6. — The abundances relative to Solar are plotted against their first ionization potential. Results from 
the current work are marked with squares, and the uncertainties are one standard deviation from the Monte 
Carlo fit. Results from low-resolution X-ray determinations are also plotted (diamonds and asterisks, with 
dotted error-bars) as well as the AR Lac photospheric values (triangles) and scaled Solar coronal abundances 
(dashed line); the symbol legend indicates the first author. 



29 




-1000 -500 500 1000 -1000 -500 500 1000 

Velocity [km/s] Velocity [km/s] 



Fig. 7. — In this figure we compare the normahzed flux profiles at quadrature and conjunction phases. 
Noniialization is done within the wavclcngtli band by subtracting tlic minimum and then dividing by the 
maximum. The sohd hne profile is from the conjunction phases when the stars' line-of-sight velocity difference 
is zero, the dashed quadrature, and the lowest dash-dot line is the difference. The left panel is from MEG 
and the right is HEG (at the same wavelength, HEG has twice the resolution of the MEG). Both broadening 
and shifts are apparent. 
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Fig. 8. — In this figure wo plot the line flux residuals normalized by their measurement uncertainty (upper 
graph) and the ratio of measured line flux to model flux (bottom graph). The abscissa is the first moment of 
the line's emissivity distribution, rather than the temperature of maximum emissivity. This is a conservative 
bias for the hydrogen-like ions, since they have a long tail on the high-temperature side they move to slightly 
higher values than their peak. In the both plots, each ion is denoted by the first character in its element 
abbreviation with some use of lower-case for duplicate characters (see the legend). While many features 
lie within ±1(7 of zero, many lie outside two and 3(T, in excess of what would be expected statistically for 
Gaussian measurement errors alone. Hence, there isMdeficiency in the DEM and abundance model, errors 
in the atomic data, mis-identified and blended lines, or likely, some combination of these. 



